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Abstract 

The Runge-Gross [E. Runge, and E. K. U. Gross, Phys. Rev. Lett., 52, 997 (1984)] 
action functional of time-dependent density-functional theory leads to a well-known 
causality paradox, i.e., a perturbation of the electronic density in the future affects 
the response of the system in the present. This paradox is known to be caused by 
an inconsistent application of the Dirac-Frenkel variational principle. In view of the 
recent solutions to this problem, the action functional employed by Runge and Gross 
in their formulation of time-dependent density functional theory is analyzed in the 
context of the Keldysh contour technique. The time-dependent electronic density, as 
well as the concept of causality, are extended to the contour. We derive a variational 
equation that obeys causality and relates the exchange-correlation potential with its 
kernel, and the functional derivative of the exchange-correlation action functional with 
respect to the density. It is shown that the adiabatic local-density approximation is a 
consistent solution of this equation and that the time-dependent optimized potential 
method can also be derived from it. The formalism presented here can be used to 
find new approximations methods to the exchange-correlation potential and avoid the 
causality dilemma. 

1 Introduction 



Time-dependent density-functional theory (TDDFT) [1-3] establishes the time- dependent 
(TD) electronic density as the primary object of study to understand the dynamics of molec- 
ular systems. TDDFT is widely used to calculate spectroscopic properties of molecules and 
solids, specially when TD perturbation theory is not applicable [4]. However, TDDFT can 
also be used to study electronic excitations in the linear regime, or predict the electronic 
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ground-state energy and density [5, 6]. The foundation of TDDFT is the theorem of Runge 
and Gross (RG) [7] stating that there is a one-to-one mapping, given an initial state, between 
electronic TD densities and TD external potentials. Later van Leeuwen [8] showed that it 
is possible to reproduce the TD electronic density of the system of interacting electrons by 
a system of non-interacting electrons, which makes possible the use of the TD Kohn-Sham 
(KS) equations. Challenges in TDDFT include the correct description of charge transfer 
excitation [9, 10], electronic transport through a molecule connected to metallic leads under 
a bias [11-13], high-order harmonic generation [14], double excitations [15], van der Waals 
interactions, among others [16]. 

In TDDFT, the prediction of the evolution of the electronic density is reformulated in 
terms of the TD KS equations, which are easier to solve than the time-dependent Schrodinger 
equation (TDSE). Moreover, every observable of the system can be expressed as a functional 
of the density because the wave function is a density-functional as well. However, Runge 
and Gross [7] proved that the one-to-one mapping between TD densities and TD potentials 
is valid under the restriction that the TD external potential is Taylor-expandable in terms 
of the time variable. The question as to how vast the set of TD potentials (or TD densities) 
of the RG theorem should be is still an open question. The proof shown by Ruggenthaler 
and van Leeuwen [17] and Ruggenthaler et al. [18] suggests that the Runge-Gross theorem 
can be extended to a wider set of TD potentials including those that are non-analytic in 
time. Despite it is known that this map exists and there is a formal procedure to construct 
KS potentials, it is still a challenge to calculate the TD potential from a given TD density 
[18-20]. 

In ground-state density-functional theory (DFT), the exchange-correlation (XC) potential 
is expressed as the functional derivative of the XC energy functional with respect to the 
time-independent electronic density. In TDDFT, an analogous variational relation between 
the TD XC potential and its action functional has been sought for the last three decades. 
Peuckert [21] first suggested that the Dirac-Frenkel action functional and its variational 
principle should be used in TDDFT. Later Runge and Gross [7] showed that the Dirac- 
Frenkel action functional extended to TDDFT leads to identify the TD XC potential as a 
functional derivative of the XC action functional with respect to the density. However, it was 
later found by Gross et al. [22] that this gives rise to a paradox in which a variation of the 
density in the future induces a perturbation of the potential in the past. According to this, 
the inverse first order response functional would not be causal. Gross et al. [22] conjectured 
that the paradox could be solved by incorporating the causality principle explicitly into the 
action formalism. 

To resolve this causality paradox several works have been published. Rajagopal [23] 
introduced an action based on the work of Jackiw and Kerman [24] in quantum field theory. 
However, this formalism does not use the density as basic variable but a transition density 
that can be negative- valued; this quantity is unsuitable as a basic quantity for TDDFT. 

van Leeuwen [25] proposed a functional in the Keldysh contour with similar properties to 
that of a free energy. This functional depends on a pseudo-density in the Keldysh space that 
reduces to the density of the system when the potential in the Keldysh space corresponds 
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to a physical potential. Due to the symmetry properties of the first-order response function 
of the pseudo-density, the causality is restored when the density is mapped to the real-time 
regime. However, the van Leeuwen formalism requires expansion of the action functional in 
terms of Feynman diagrams, while the functional of RG does not require such expansion. 
Furthermore, the operator used in this formalism for the pseudo- density is not Hermitian in 
general and thus the pseudo-density does not integrate to the total number of electrons of 
the system; except when the density is physical. 

Recently, Vignale [26] solved the causality paradox in real-time by showing that the 
source of the problem in the RG formulation is a boundary condition. He showed that 
only the initial condition is necessary in the Runge-Gross functional to recover the causality 
restriction in general, and derived an expression for the XC potential that is causal. 

In this paper I review Vignale's solution of the causality paradox in real-time from the 
perspective of unitary propagation and later use Vignale's theory to extend the RG action 
functional to the Keldysh contour. The RG action functional in the Keldysh space, unlike 
the van Leeuwen functional, does not require diagrammatic expansion and uses an electronic 
density that is a causal functional of the potential in the Keldysh space. By the RG theorem 
applied to the Keldysh space and under the assumption that the density is a strictly causal 
functional of the potential in the contour, I show that a variational equation relating the XC 
potential with the XC action functional arises. This equation shows an explicit dependence 
on the memory of the system through the XC kernel. I show that the adiabatic local density 
approximation (ALDA) is consistent with this equation, how the TD optimized potential 
method (TDOPM) can be obtained, and also how the ground-state XC potential can be 
recovered. 

2 Real-time Analysis 

The Dirac-Frenkel variational principle provides a method to derive the TDSE and its ap- 
proximations by finding a stationary value of the action functional: 



where T and W are the kinetic energy and electron-electron repulsion energy operators, 
respectively, and n(r) is the density operator. The Dirac-Frenkel functional is defined over 
a Hilbert space of antisymmetric wave functions representing bound systems of N electrons. 
The TD Schrodinger equation (TDSE) is thus obtained by setting 




(1) 



In this work we consider Hamiltonians of the form: 




(2) 



SfWlfav] = . 



(3) 
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This variational principle supposes that 5ip(t ) = 8ijj{tx) = 0. The solution of this equation, 
denoted as ^[v](i), is the solution of the TDSE: 

i^\m(t)) = H[v)m[v)(t)) ■ w 

is said to be a f-representable ket in real-time, which expressed in terms of the 
unitary evolution operator is: 

\iP[v](t)) = U[v)(t,to)m )} , (5) 

where t 

U[v](t,t ) = texp(-i£ ds #[«](*)) . (6) 

Here T is the time-ordering operator in real-time. In this work we interpret the integral in 
the above equation to be taken over the interval \po,t), i.e. 

ft rt-\A 

/ ds H[v)(s)ds = lim / ds H[v)(s) . (7) 

J to t J to 

Note that the ket IV'M^)) is a causal functional of the potential: It is determined by the 
potential v at times less than t. We refer to this dependency on the potential (Eq (6)) as 
the strict causality assumption. 

All observables of ip are also causal functionals of u. For example, the density of the 
system, 

n[v](rt) = ^(to)|^MMo)M^M(Mo)|^o)> , (8) 

is determined by v in the interval [to,t) [3]. By the RG theorem, given a fixed initial state, 
the potential v at times in [to,t) uniquely determines n in the interval [to,t), and vice versa. 
If we denote as u[n] the external potential as a functional of the TD density, then a first 
order variation in u is given by a variation of n over the interval [to, t): 

5u[n](rt)= [ dt' / dV x^Hirt^r't^nir't') , (9) 

J to J 

where 

This indicates that 5u{vt) / 5n{r't') for t < t' is not defined because it does not contribute to 
the integral of Eq. (9). However, for convenience we set: 

X" 1 [n](rt,r , t') = * < ^ • ( n ) 

Vignale [26], however, employing the evolution equation of the current, showed that 
5u(rt) /5n(r't r ) is related to 5(t — t') and its first and second order time-derivatives when 
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t = t'. This result is obtained under two assumptions different from ours: First, u(rt) is 
determined by n(r't') for t' < t. And secondly, the functional derivative of the stress tensor 
with respect to the density vanishes at equal times. Our assumption avoids this singularity 
in x -1 and will be used to simplify our calculations in the Keldysh space (section 3). 
Now let us consider the Runge-Gross action functional: 

A v [n] = ^ dt (VM(i)|i| - H[v](t)mn](t)) , (12) 

where |-?/>[n](t)) = [n]] (t)) , v is some TD external potential, and t 1 > t . Note that the 
ket IV'W(O) ^ s causai ; i- e -> it is determined by n in the interval [t ,t). Runge and Gross [7], 
based on the Dirac-Frenkel variational principle, imposed 8ip(to) = 5ip(t\) = and 

<U " , (13) 



5n(rt) 

which leads to the following alternative form of the variational principle 

5 B\n 



5n(rt) 

where B[n] is the internal action: 



v(rt) = , (14) 



B[n] = I dt (VN(t)|i^ - f - W\j>[n]{t)) . (15) 



tQ 

If Eq. (14) were valid then we could assert that [27]: 

r , 5 Bin] , , 

= ^ (16) 

Unfortunately, when the above function is further differentiate with respect to n, 

. ,. 5 2 B\n] 

*< rt - ri > = M rt -)L( ri ) ' (17) 

one finds an inconsistency because the above equation implies that x( r ^ r '0 7^ for t < t' . 
This is known as the causality paradox [22]. The solution to the paradox was found by 
Vignale [26], who pointed out that, according to the definition of f-representable wave- 
function, we can only set 5ip(to) = because a perturbation 5n(rt), in general, will induce a 
response 5ij)(t\) ^ 0. 

The solution of Vignale [26] can be viewed as a direct implementation of the causality 
principle into the RG functional. For example, the internal action B[n], using the TDSE, 
can be written as [3]: 



B[n} = jf" dt J d 3 r u[n\{xt)n{rt) . 
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The density-functional u[n] is causal by the RG theorem. If we differentiate the above 
functional with respect to the density and insert the result into the functional derivative of 
the RG action functional we obtain 



5A V 



u[n\{vt) - v(rt) + J dt' j dV x _1 N(rt,r't>(rY) . (19) 



5n(rt) 

Now let n v be the TD density corresponding to v , then: 



5 A: 



8n(rt) 



dt' J dV x~ l [nv\(?t,Y't')n v (Y't') . (20) 



This last equation is an alternative form of the Vignale variational formulation that shows 
that n v is not a stationary value of A,[n]. This is a consequence of constraining the wave- 
functions of the RG functional to be density-functionals of the form ^[it[n]]. Ruggenthaler 
and Leeuwen [27] showed that not every TD wave-function can be associated with a TD 
external potential (or a TD density). Hence the domain of the RG functional is just a subset 
of the domain of the Dirac-Frenkel functional, explaining why the RG and the Dirac-Frenkel 
functionals lead to different results. 



3 Keldysh-space Analysis 

Let us add a super index + or — to the time variable t. The Keldysh contour, C, is expressed 
as C = C + U C~ , where C + = [to,tf] and C~ = [t^ , We denote z as a variable in the 
contour C, and let z- x = t$ and Zf — t^. The arrow of time in C points from to and 
from ti to t$. Thus, any z G C~ is said to be later than any z' G C + . If z, z' G C~ we say 
that z is later than z' if t(z) < t(z'), where t(z) is the real value of z. A ket in C is denoted 
as |^ c [u c ](^)), where u c (rz) is some potential in C. A physical potential in C is denoted as u c 
and it satisfies u c (rt + ) = u c (rt~). Thus a potential in real-time is mapped to the Keldysh 
space when setting ■u c (rt ± ) = u(rt) (t^ we denotes evaluation at C + or C~). 
We now extend the unitary propagator U to the Keldysh space as follows: 



U c [u c ](z,Zi) = Tcexp - ij dz' H c [u c ](z') , (21) 



where 7c is the path-ordering operator in C (for example, l~c[B c {z')A c (z)\ = A c (z)B c (z') if 
z is later than z'). The Hamiltonian in the Keldysh space now reads H c [u c ](z) = T + W + 
J d 3 r u c (rz)n(r). The integration over the pseudo-time is defined as: 

u ' ~ U: + /«, m un, * = r. (22) 

A w-representable ket in C is thus expressed as |-?/' c [m c ](z)) = U c [u c ](z, Zi)\ip c (zi)) , where 
\^c{ z i)) = l^(*o)) is the initial state of the system. Note that ip c (z) does not depend on the 
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potential u c at later times than z. As in the real-time case, we assume that the end point of 
the integral in Eq. (21) is not included; this can be considered as an extension of the strict 
causality assumption to the Keldysh contour. We define the density in C as [28]: 

n c [u c ](rz) = (Ul[u c )(z, Zi)h(r)U[u c ) c (z, Zi)) , (23) 

where (•) = (ip c (z[)\ ■ \ip c (z\)). To prove that there is a one-to-one mapping between n c and 
u c , it is sufficient to notice that ip c satisfies the Schrodinger equation in C + . Therefore, if 
the potential can be expressed as a power series around z h then the RG theorem and its 
extension [17] including non-analytic potentials apply in this case. 

Let us examine the action functional proposed by van Leeuwen [25], which reads 

A vL [u c } = iln(U c [u c }{ Zi , Zi )) . (24) 

The functional derivative of this functional with respect to the potential u c yields the pseudo- 
density [25]: 

n vL (rz) = (UMl^m^l . (25 ) 
(U c (z h Zi)) 

However, the above density is an average of the operator: 

^vL,H(r) = U c (z{, z)n(r)U c (z, Zi) , (26) 
which is not a Hermitian operator. Therefore 

d 3 r n vLiH (vz) = NU c (z uZi ) , (27) 



where N is the particle-number operator. This implies that n v L does not integrate to N\ 
except when the potential u c is physical [25]. The density n c , on the other hand, integrates 
to iV and is always positive. 

It can be shown that the response function of the density in C is given by: 

Xc [u c ](rz,r'z') = 6u c pjj = -i([r^ CI HK](r^),n c ,HK](r / 2; , )]) , (28) 

where the Heisenberg representation of the density operator n(r) is 

n c ,n[u c ](rz) = Ul[u c ](z, Zi)n(r)U c [u c ](z, Zi) . (29) 

Eq. (28) is valid if z is later than z', and we set Xc{ rz , r ' z ') — if z' is later than or equal 
to z. 

The inverse first order response function x~ 1 [ra c ](r2;, r'z') = 5u c (rz)/5n c (r'z'), according 
to the RG theorem extended to the C, must also satisfy causality in the contour, e.g., 
X c T 1 (r2;, rz') = if z = z' or z' is later than z. When a physical potential is used, the 
Heinsenberg operators recover their usual form in real-time. Therefore, we obtain a physical 
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density 77, c (rt ± ) = n(rt). From Eq. (28) we can show that the first order response function 
satisfies the antisymmetry relationship: 



Xc (Tt + ,r't' + 



-Xc(r't'-,rr 



(30) 



where u c = u c denotes evaluation at the physical regime. Note that Xc also satisfies 
Xc(rt+, r't'+) = Xc(rt~, r'f+) and Xc(r'f~ , rt + ) = Xc(r'f~ , rt~) if t > t' and u c = u c . 
The response of the density in the Keldysh space is [25]: 



5n c [u c ]{rz) = J dz' J d 3 r' Xc[u c ){rz, r'z')6u c (r'z' 



(31) 



To obtain the response in real-time, the variation of a physical potential must satisfy 
Su c (rt + ) = 5u c (rt~) = Su(rt). Using the aforementioned properties of \c to calculate the 
above integral, the response of the density turns out to be independent of the time location 
in the contour, i.e. 5n c (rt + ) = 5n c (rt~) = 5n(rt). Hence, it is determined by: 



5n c [u c }(rt) = jf" dt' J dV X c[u c }(vt ± ,v't' + )5u c (v't l+ ) 



This result allows us to identify the response in real-time x( r ^ r 't') as Xc( r ^ ± ) r 't H 
—Xcijc't'^, rt~)\ Uc=ilc , which are causal. Exchanging variables in the integral of XcXc 1 reveals 
that Xc 1 satisfies the same relationships of Xc regarding exchange of variables at physical 
densities. 

Let us extend the functional A v to the Keldysh space: 



(32) 



or 



Av c [n c ] = B[n c ] - J dz J d 3 r n c (rz)v c (rz) 



where 



B\n c 



H Q 

dz (V'cKK^Ii^- 



f-WMn c )(z)) , 



(33) 



(34) 



v c is some external physical potential, and df(z)/dz = df(t a )/dt, where a = +, — . Vignale 
equation in this case reads: 



513 



§n c (rz) 



VAYZ) 



■I i ( M^cfek 



Sn c (rz) 



J ' ^' J d 3 r' n c ,, c (rV) Xc - 1 K, c ](r'z',rz) 



(35) 



The left hand side of the above equation corresponds to 8Av c /Sn c (rz) evaluated at the density 
that yields v c , n c ,v c - Additionally, the above equation also gives the functional derivative 
5B/5n c (rz) for an arbitrary density n c ; in this case, we replace v c by u c [n c ](rz), n c ^ c by n c , 
and the inverse response function has to be evaluated at n c . 
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Let us introduce the KS action functional: 

-Zf 



A,s c>s K] = B s [n c ) - J dz J d 3 r n c {rz)v C)S {rz) , (36) 
where v CjS (rz) is some effective external potential and 

B s [n c ] = £ dz ($ C)S [n c ](z)|i^ - f |$ c > c ](z)) . (37) 
The KS wave function is a Slater determinant of TD KS orbitals {<fi c ^(rz)} that satisfy: 

= ( - + u c>B [n c ](vz)y c!i (vz) , (38) 

where m c , s [^c] is the KS potential that represents n c {rz). Thus, if we differentiate B s with 
respect to the TD density we obtain: 



5n r (rz) 



u c , s [n c ](rz) + J dz' J d 3 r' n c (r' z')x c ^[n c \{r' z' , rz) 



(39) 



where x c I ( YZ i r ' z ') — Su C) g(T r z')/Sn c (Tz). 
Recall the Hartree functional: 



, r n 1 T , f ,q f , q , n c (r'z)n c (rz) , , 

^H[n c ] = - / dz dh dV c \ ; c , ; . (40) 



2 . 

Let us introduce the XC action functional: 

Ache] = B s [n c ] - B[n c ] - A H [n c ] . (41) 

Using Eqs. (35) and (39) it is found that the functional derivative of the XC action functional 
can be expressed as: 

« c ,xcM+ / dz' f dh'n c (r'z')[x~l(r'z',rz) 

6A (42) 

-X c (r^,r2;)J 



Here M CjXC (rz) = u c ^{rz) — u c (rz) — u C} k(tz), where the Hartree potential is u Ci n[n c ](rz) = 
J d 3 r'n c (r'z)/|r — r'|. Now introduce the XC kernel f C)XC (rz,r' z') = Su C;XC (rz)/Sn c (r'z'), 
which satisfies: 

X:}(rz, r'z') = x:\rz, r'z') + + f^TZ, r'z') . (43) 

The delta function in C space is defined such that JJ dz' f c (z')S c (z — z') = f c (z). The KS 
response function and the XC kernel satisfy the same properties of Xc regarding exchange of 
variables. 
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In order to simplify Eq. (42), suppose that the density is physical, n c = n c . This imposes 
that the XC potential is the same in both C + and C~ spaces. For example, if z = t + then 
the integral in time can be split up into two integrals: The first one runs from t + to t~ , and 
the second one from t~ to £q . There is no contribution from the first integral due to the 
symmetry properties of Xc 1 and Xcs at physical densities. For the second integral we can 
use the antisymmetry relation to obtain in real-time that: 



lly Jrt)+ I df l dV U(rt, r't')n(r't') = - ( 44 ) 



(rt)+J t dt'J 
where u xc (rt) = ■u CiXC (rt ± ) and 

SA XC 5A 



dn(rt) 



dn(rt) 5n c (rt±) 



(45) 



Setting z = t~ in Eq. (42) also leads to Eq. (44) when n c = n c ; for this reason we 
expressed the final result in real-time. Because / CiXC in the C space also has the same prop- 
erties as Xc 1 we identify the XC kernel in real-time, f xc (rt,r't'), as / c , X c( r * ± , r't /+ )\ nc =n c , or 
— / c ,xc(r'£' - , rt~)\ nc=nc . Thus, the XC kernel is causal in real-time. 

Given that we assumed that the response functions Xc and Xc,s are strictly causal in C, 
the integral in Eq. (42) is taken over the interval (z, zt\. This implies that the Hartree kernel 
5 c (z — z')/\r — r'| lies outside the integration limits and thus it has no contribution to Eq. 
(44). Based on this, the integral in Eq. (44) is carried out strictly over the past of t, i.e., 
[to,t). Hence, our causality assumption avoids singularities at equal-times and simplifies the 
transition to real-time. 

Eq. (44) is the main result of this work. It is a variational equation that establishes 
a causal connection between u xc in real-time with an XC action functional in the Keldysh 
space, and the memory of the system. If an approximation to the XC action functional is 
known, then Eq. (44) can be used to estimate the XC potential. The potentials u(rt) and 
u s (rt) also satisfy the same type of equation as that of u xc ; one has to replace / xc and A xc 
by X~ l and £>, or Xs 1 and B s . 

Note that the left-hand side of Eq. (44) is a functional of the density n c (rt ± ), or simply 
n(rt). This implies that the second functional derivative of A xc with respect to the density 
in real-time is not symmetric, i.e.: 

- SAxc =o t' > t (46) 

Sn(r>t>) Sn(rt) ~ ' 1 ' 

Here, the symbol 5/5n(r't') represents regular functional differentiation in real-time because 
the operation S/Sn(rt) already involves evaluation at the physical regime. The above result 
is a consequence of implementing causality in the C space explicitly using the path-ordering 
operator. Furthermore, recursive differentiation of Eq. (44) also allows us to express its 
solution as a series of functional derivatives of A xc . This reads 

dA xc 

«xc(xi) = ir~r~\ + w xc(xi) (47) 

on(xi) 
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where 

Wxc(xi) = ^ 2 J d ^ X2) ' ' ' dMxm W m )---Mx 2 )Hxi) • (48) 

Here x m = r m , t m , m = 1, 2, . . . and d/i(x m ) = n(x m )d 4 x m . The functional derivatives in the 
integral are zero if, for any i > j,U>tj. This series shows that the XC potential depends on 
perturbations of the XC potentials in all orders. However, in order to achieve convergence 
the functional derivatives must decrease as their order increases. 

Now let us apply our variational equation to the derivation of the ALDA XC potential: 



A^ K [n c } = J\z j d 3 r [e xc (n)n] 



(49) 

n=n c (rz) 



where e xc is the local XC energy density. To solve Eq. (44) the memory term can be neglected 
to yield 



(50) 

n=n c (rz) 



Further differentiation leads to the kernel formula: 

f ALDA ' - •- ~> '- a 



f: c WA (rt, r't') = 5(r - r')5(t - t')— 2 [e xc (n)n] 



(51) 

n=n(rt) 



The singularity of the XC kernel does not contribute to the integral term of Eq. (44) because 
the end point is not included, or in other words, the end point is approached in a limiting 
procedure. Hence, the above equation satisfies Eq. (44) and thus it is the solution of 
it. The singularity of the XC kernel arises from the definition of the XC potential, which 
implies that at equal-times the XC kernel must cancel the singularity of the Hartree kernel. 
However, the ALDA XC kernel does not cancel the singularity of the Hartree kernel due to 
the self-interaction error. 

Another application is the TDOPM. The exchange functional form remains the same as 
the one proposed by van Leeuwen [25]: 

4k] = £ dz (*c[n c ](z)\W\*c[n c ](z)) - An[n c ] . (52) 

To derive the TDOPM one has to assume that (for example, see [3]): 

5Av r t^W.s 



8n{rz) 5n(rz) 



(53) 



If we set A xc = A x and expand A^ c using Eq. (41) we find that the memory term in Eq. 
(44) can be discarded. Hence we can write: 
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The right hand side of the above equation can be calculated using the chain rule. If the 
result is multiplied by Xs and then integrated, the final result coincides with that of Ullrich 
et al. [29]. 

Ground-state DFT is also accessible with this theory. We can introduce a slowly varying 
density nj(rz) = n c (rz/T), where T — > oo. One can use the adiabatic theorem to show that: 

rzf 

lim A xc [n T c ] = lim / dz E xc [n T c (■ , z)} , (55) 

T— >oo T->-oo J 

where E xc is the XC energy functional of DFT. The above equation is local in time. As in 
the previous case, the solution of Eq. (44) has to be of the form: 

(56) 

n=n T (rt) 

where n T = n^. 

4 Conclusions 

To summarize, we examined the RG action functional and the solution of the causality para- 
dox by Vignale [26] from the point of view of unitary evolution. We extended this solution 
to the Keldysh space, and, under the strict causality assumption, we found a variational 
equation for the XC potential that involves an XC memory term. The solution of this vari- 
ational equation is a series in terms of functional derivatives of the XC action functional 
in the Keldysh space. We showed that it is possible to derive the ALDA XC and TDOPM 
exchange potentials from the present theory and that ground states are also accessible using 
the adiabatic theorem. 
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